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LONG-TERM  GOALS 

Modeling  and  predicting  ocean  optical  properties  for  coastal  waters  requires  linking  optical  properties 
with  the  physical,  chemical,  and  biological  processes  in  the  upper  ocean.  Our  long-term  goal  is  to 
incorporate  optical  processes  into  coupled  physical-biological  models  for  coastal  waters,  develop  and 
improve  integrated  ocean  forecasting  systems,  including  prediction  of  ocean  optical  properties. 

OBJECTIVES 

1)  To  improve  performance  of  the  coupled  physical-biological  model,  which  is  based  on  the  Regional 
Ocean  Model  System  (ROMS),  for  the  Monterey  Bay; 

2)  To  incorporate  optical  variables  into  a  one-dimensional  (ID)  physical-biological  as  well  as  the 
improved  coupled  3D  physical-biological  model  for  the  Monterey  Bay; 

3)  To  use  these  variable  to  drive  a  radiative  transfer  model  (HydroLight/EcoLight)  that  will  simulate 
and  predict  the  subsurface  light  field  as  well  as  the  ocean’s  color; 

4)  To  conduct  ROMS-biological-optical  model  simulations  for  the  Monterey  Bay  with  focus  on  two 
intensive  field  experiments,  ANOS  II  (August  2003)  and  Monterey  Bay  2006  (August  2006),  and  use 
the  bio-optical  measurements  to  constrain  and  improve  the  models. 
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APPROACH 


Our  approach  is  to  incorporate  spectrally-resolved  inherent  and  apparent  optical  properties  into  a  one¬ 
dimensional  (ID)  multi-nutrient,  phytoplankton,  zooplankton  and  detritus  ecosystem  model.  By 
comparing  the  bio-optical  model  output  with  in  situ  and  satellite  ocean  color  measurements,  we 
highlight  the  fact  that  optical  properties  play  an  important  role  in  identifying  and  reducing  uncertainties 
in  ecosystem  models,  which  provide  constraints  for  determining  variables  and  related  parameters.  Our 
long  term  goal  is  to  incorporate  the  improved  ID  bio-optical  modeling  into  the  Regional  Ocean  Model 
System  (ROMS)  for  the  Monterey  Bay,  and  make  prediction  of  coastal  ocean’s  optical  properties  on 
time  scales  of  1  to  5  days,  i.e.,  the  time  scales  of  accurate  atmospheric  forecasts. 

The  modeling  system  constructed  in  this  study  consists  of  four  individual  models,  physical-ecosystem 
model,  photo-acclimation  model,  optical  model,  and  radiative  transfer  model.  The  physical-ecosystem 
model  used  in  this  study  is  based  on  the  Carbon,  Silicon,  Nitrogen  Ecosystem  (CoSINE)  model  (Chai 
et  al.,  2002),  which  was  originally  developed  to  simulate  biogeochemistry  in  the  equatorial  Pacific 
upwelling  region.  To  reproduce  observed  variation  in  phytoplankton  chlorophyll  to  carbon  ratio  with 
growth  conditions  (light,  nutrients  and  temperature),  Geider  et  al.  (1987,  1998)  developed  a  photo¬ 
acclimation  model  with  single  nutrient  of  nitrogen.  Moore  et  al.  (2002)  modified  the  photo-acclimation 
model  so  that  the  model  could  be  embedded  in  multi-nutrient  ecosystem  model,  which  was  used  in  this 
study.  We  developed  an  optical  model  that  explicitly  represents  spectrally-resolved  IOPs  (absorption, 
scattering  and  attenuation),  apparent  optical  properties  (AOPs,  such  as  diffuse  attenuation),  and 
radiometric  quantities  such  as  PAR.  The  optical  model  requires  the  inputs  from  the  physical-ecosystem 
model,  such  as  phytoplankton  associated  chlorophyll  and  non-algal  particles  (NAP).  The  radiative 
transfer  model  (Ecolight;  Sequoia  Scientific,  Inc.)  provides  downwelling  irradiance  at  sea  surface  and 
underwater  average  cosines  of  downwelling  light  (from  400  to  700  nm),  and  the  hourly  values  are  used 
as  external  forcing  for  the  optical  model.  The  physical-ecosystem  model,  photo-acclimation  model, 
optical  model  and  radiative  transfer  model  are  connected  as  one  bio-optical  modeling  system.  We 
applied  this  system  to  the  equatorial  Pacific  upwelling  region  (5°S-5°N,  90°-180°W,  the  “Wyrtki 
Box”),  where  the  physical-ecosystem  model  originally  developed  (Chai  et  al.,  2002)  and  some  optical 
properties  have  been  measured  for  this  region  (Gardner  et  al.,  2003;  Behrenfeld  et  al.,  2006). 

WORK  COMPLETED 

This  report  summarizes  modeling  activities  between  1  October  2006  and  30  September  2007. 

We  have  completed  the  equatorial  modeling  effort.  This  includes: 

1 .  Developed  a  bio-optical  module  that  translate  the  biogeochemistry  into  inherent  optical 
properties  (IOPs). 

2.  Incorporated  a  radiative-transfer  model  that  predicts  the  subsurface  light  field  given  the  IOPs. 

3.  Performed  a  sensitivity  study  to  the  parameters  of  the  bio-optical  module  based  on  their 
uncertainties  (see  Table  1  below). 

4.  Compared  our  results  with  IOPs  measured  in  situ. 

5.  Presented  our  results  in  the  Ocean  Optics  conference  and  had  our  results  peer-reviewed  and 
published  in  the  journal  Biogeosciences  (Fujii  et  al.,  2007). 

6.  Assembled  a  dataset  of  all  the  IOPs  in  Monterey  Bay,  which  will  be  the  site  of  our  next  modeling 
effort. 
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RESULTS 


Given  that  in  our  past  reports  we  focused  on  the  matchup  between  modeled  and  observed  properties 
here  we  focus  on  a  sensitivity  analysis  to  the  bio-optical  model.  Such  analysis,  to  our  knowledge,  has 
never  been  perfonned  and  provided  many  lessons  for  us  regarding  what  is  and  what  is  not  essential  to 
include  in  the  model. 

The  values  of  parameter  chosen  for  this  study  are  based  on  observations  but  the  observed  values  have 
substantial  variability  that  arises  from  environmental  and  methodological  variability.  To  elucidate  how 
model  results  are  affected  by  variations  in  the  optical  model  parameters,  we  conducted  a  sensitivity 
study  of  the  model  to  those  parameters  by  changing  their  values  individually  by  30%  of  their  standard 
values.  Such  variability  encompasses  the  bulk  of  observed  values  (Table  1). 


Table  1.  Sensitivity  of  model  results  to  optical  parameters:  Carbon-specific  absorption  coefficient  by 
NAP  at  44 Onm  (a+ nap(440)),  absorption  coefficient  by  CDOM  at  440nm  (uCdom(440)),  background 
backscattering  coefficient  (bhbg),  ratio  of  phytoplankton  carbon  to  POC  (Rv poc),  backscattering  ratio 
for  picoplankton  (bh pi),  diatoms  (bh _P2),  NAP  (bh nap),  and  background  particles  ( bhbbg )•  Values  in 
parentheses  denote  model  results  from  a  control  run  with  bio-optical  model.  Euphotic  layer  depth  is 
defined  as  a  depth  of  0.1  %  of  sea  surface.  Sources  noted  here  are:  (1)  Babin  et  al.,  2003a  and 
2003b;  (2)  Simeon  et  al.,  2003;  (3)  Stramski  and  Kiefer,  1991,  Clio  and  Azam,  1990,  and  Behrenfeld 
and  Boss,  2006;  (4)  Eppley  et  al.,  1992,  DuRand  et  al.,  2001,  Gundersen  et  al.,  2001,  and  Oubelkheir 
et  al.,  2005;  (5)  Twardowski  et  al.,  2001;  (6)  Twardowski  et  al.,  2001,  and  Boss  et  al.,  2004. 
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We  found  the  model  results  of  biogeochemical  properties,  i.e.,  surface  NO3,  Si(OH)4,  and  maximum 
chlorophyll  and  its  associated  depth,  to  be  most  sensitive  to  changes  in  the  absorption  coefficient  by 
CDOM  at  440nm  (acDOM(440)).  The  surface  NO3  increases  and  Si(OH)4  decreases  with  the  increase  of 
acDOM(440),  due  to  an  increase  in  contribution  by  diatoms  when  acDOM(440)  is  higher.  The  maximum 
chlorophyll  decreases  and  appears  at  a  deeper  layer  of  85m  with  an  acDOM(440)  increase.  These  model 
results  reveal  that  the  CDOM  concentration  strongly  affects  phytoplankton  community  structure  and  its 
dynamics.  While  CDOM’s  inherent  effects  on  backscattering  and  hence  beam  attenuation  coefficient  at 
660nm  (Cp(660))  are  negligible,  an  increase  of  acDOM(440)  yields  a  Cp(660)  decrease  due  to  a  decrease 
in  small  algal  POC.  The  modeled  euphotic  layer  depth,  defined  as  a  depth  of  0.1%  light  level  of  sea 
surface,  decreases  from  120m  to  1 10m  by  varying  acDOM(440)  from  0.01 1  to  0.021  (m"1),  primarily  as 
a  result  of  enhanced  absorption  by  CDOM.  The  change  of  the  euphotic  layer  depth  is  relatively  small 
because  PAR  is  more  controlled  by  absorption  by  water  than  absorption  by  underwater  particle  and 
CDOM  concentration.  However,  the  euphotic  layer  depth  is  more  sensitive  to  acDOM(440)  than  to  the 
other  optical  parameters  due  to  significant  absorption  by  CDOM  at  short  wavelengths  around  400nm, 
at  which  absorption  by  water  is  negligible. 

The  observed  sensitivity  to  CDOM  concentration  is  the  result  of  CDOM  absorbing  light  that  would 
otherwise  be  absorbed  by  phytoplankton.  This  effect  is  more  pronounced  on  picoplankton  as  they  have 
a  relatively  higher  portion  of  their  energy  absorbed  in  the  blue  wavelength  where  CDOM  absorbs  (they 
are  less  packaged  and  thus  have  a  higher  blue  to  red  absorption  ratio).  Changes  in  the  relative 
abundance  of  small  and  large  phytoplankton  results  in  a  change  in  the  biogeochemical  properties  of  the 
upper  ocean  since  their  metabolic  requirements  and  interaction  with  other  trophic  levels  are  different. 
Variation  in  other  optical  parameters  also  contributed  to  changes  in  the  model  results,  but  their  impact 
is  smaller  than  that  of  acDOM(440).  The  carbon-specific  absorption  coefficient  by  NAP  at  440nm 
(a+NAp(440))  has  weaker  but  similar  effects  on  surface  nutrient  and  maximum  chlorophyll 
concentration  as  acDOM(440).  The  modeled  aNAp/ap(440)  is  the  most  sensitive  to  a+NAp(440),  changing 
by  a  factor  of  2  from  0. 10  to  0.20.  The  modeled  Cp(660)  is  affected  most  significantly  by  the  ratio  of 
phytoplankton  to  particulate  organic  carbon  (R'Vx).  The  modeled  euphotic  layer  depth  does  not 
change  from  the  standard  value  of  1 15m  by  changing  any  optical  parameters  except  for  a+NAp(440), 
which  indicates  the  important  role  of  absorption  in  the  lower  layer  below  chlorophyll  maximum  in 
determining  the  euphotic  layer  depth.  Varying  any  of  the  backscattering  ratios,  regardless  of  particle 
type,  does  not  affect  the  modeled  biogeochemical  properties  but  does  influence  Cp(660).  However,  the 
sensitivity  of  Cp(660)  to  the  backscattering  ratios  depends  on  particle  type,  being  stronger  for  NAP  and 
picoplankton  and  weaker  for  diatoms  and  background  particles,  reflecting  the  higher  backscattering 
coefficient  by  picoplankton  than  by  diatoms. 

The  overall  sensitivity  study  shows  that  narrowing  the  observed  ranges  of  optical  parameters  above  is 
required  to  reduce  uncertainties  in  reproducing  biogeochemical  properties.  In  addition,  the  above 
analysis  suggests  that  although  the  dynamics  of  neither  CDOM  nor  bacteria  are  currently  incorporated 
explicitly  in  the  model,  embedding  CDOM  as  a  state  variable  in  the  ecosystem  model  should  be  given 
priority  over  bacteria  to  improve  simulating  bio-optical  interactions.  We  will  use  these  results  in  our 
future  development  of  the  model. 
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We  developed  an  ecosystem  model  that  explicitly  represents  biogeochemically  and  optically  two 
phytoplankton  and  two  zooplankton  functional  groups,  as  well  as  multiple  nutrients  and  non-algal 
particles  (NAP).  We  applied  the  model  to  the  equatorial  Pacific  upwelling  region  and  found  that 
utilizing  an  optical  model  to  convert  from  ecosystem  model  state  variables  to  optical  parameters  and  a 
realistic  subsurface  light  provides:  (1)  more  data  to  compare  model  output  with  providing  a  more 
rigorous  test  on  model  formulation  and  choice  of  parameter  values,  especially  for  those  that  are 
difficult  to  measure  in  high  resolution  in  time  and  space,  (2)  the  required  input  to  obtain  a  realistic 
subsurface  light  field  by  linking  the  optics  to  a  radiative-transfer  model  (Ecolight),  and  (3)  improved 
simulation  realism  with  respect  to  key  biogeochemical  processes,  such  as  photosynthesis,  which  are 
crucial  for  assessing  oceanic  carbon  cycling  and  food  web  dynamics.  The  additional  optical 
measurements,  being  routinely  available  from  research  vessels,  autonomous  platforms,  and  space- 
home  observations,  can  now  be  used  directly  for  comparison  and  testing  of  the  output  of  our  new 
coupled  bio-optical  model.  This  is  an  improvement  over  the  limited  number  of  variables  that  can  be 
used  to  test  our  previous  ecosystem  models  with  no  explicit  optical  properties.  Model  sensitivity 
studies  on  optical  parameters  suggest  that  CDOM  may  have  an  important  role  in  phytoplankton 
dynamics,  nutrient  cycling,  and  light  field  in  the  euphotic  layer. 

In  preparation  of  transitioning  to  a  3-D  modeling  effort  of  Monterey  Bay  we  have  assembled  a  large 
data  set  of  biochemical  and  bio-optical  data  (for  spatial  coverage,  see  Fig.  2).  This  is  one  of  the  densest 
bio-optical  data  ever  collected  that  will  be  used  to  refine  our  bio-optical  model  and  test  our  results 
against.  We  plan  to  collaborate  with  the  scientists  involved  in  collecting  these  data  for  modeling  and 
data  analysis  for  the  Monterey  Bay  region. 

IMPACT/APPLICATIONS 

Incorporating  ocean  optical  processes  into  coupled  physical-biological  models  will  enable  us  to 
simulate  and  forecast  optical  properties  in  coastal  waters.  With  demonstration  of  some  initial  successes 
of  developing  physical-biological-optical  modeling  and  data  assimilation  capability  for  Monterey  Bay, 
we  should  be  able  to  start  the  development  of  an  end-to-end  ocean  forecasting  system.  Such  modeling 
system  would  be  a  powerful  tool  to  design  the  adaptive  sampling  strategy  and  would  be  an  essential 
component  of  future  field  experiments  both  in  and  outside  Monterey  Bay. 

TRANSITIONS 

We  have  worked  with  Drs.  John  Kindle,  Igor  Shulman,  and  Bred  Penta  at  the  Naval  Research 
Laboratory  (NRL).  The  ecosystem  model  code  has  been  transferred  to  the  Dynamics  of  Coupled 
Processes  group  led  by  John  Kindle  at  the  NRL.  The  ecosystem  model  has  been  implemented  into  both 
Navy  Coastal  Ocean  Model  (NCOM)  and  Princeton  Ocean  Model  (POM)  for  the  west  cost  of  U.S.  and 
Monterey  Bay  region.  Fei  Chai  traveled  to  the  NRL  in  late  August  of  2005  to  work  on  some  issues 
regarding  to  the  ecosystem  component  in  the  NCOM.  Fei  Chai  will  visit  NRL  in  November  2007  for 
the  ongoing  and  further  collaboration  between  the  University  of  Maine  and  NRL.  We  have  regular 
communications  between  the  NRL  and  the  University  of  Maine  on  this  collaboration.  The  improved 
ecosystem  model  and  the  optical  module  will  be  transited  to  the  NRL  Dynamics  of  Coupled  Processes 
group. 
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RELATED  PROJECTS 


This  project  has  strong  collaboration  with  other  ONR  supported  projects.  Besides  working  closely  with 

the  modeling  group  at  the  NRL,  we  are  collaborating  with  Dr.  C.  Mobley  of  Sequoia  Scientific  Inc.  for 

better  linking  the  optical  module  within  the  ROMS.  We  are  also  collaborating  with  scientists  at  the 

Monterey  Bay  Aquarium  Research  Institute  (MBARI)  for  using  the  observational  data  for  the  region. 

Yi  Chao  is  one  of  the  investigators  involved  in  the  Monterey  Bay  2003  and  2006  Experiments. 
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Figure  1:  Cruise  tracks  and  platforms  that  collected  biological  and  optical 
information  for  the  Monterey  Bay.  Optical  data  include:  chlorophyll 
fluorescence,  fluorescence  at  676nm,  optical  backscatter  at  676nm,  beam 
transmission,  beam  attenuation  coefficient,  absorption  coefficients,  PAR,  and  many  more. 
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